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Abstract. We give an algorithm that finds a sequence of approximations with 
Dirichlet coefficients bounded by a constant only depending on the dimension. 
The algorithm uses the LLL-algorithm for lattice basis reduction. We present 
a version of the algorithm that runs in polynomial time of the input. 



1. Introduction 

The regular continued fraction algorithm is a classical algorithm to find good ap- 
proximations by rationals for a given number a £ R\Q. Dirichlet [13) proved that 
for every a e M \ Q there are infinitely many integers q such that 

(1) \\qa\\<q-\ 

where ||a;|| denotes the distance between x and the nearest integer. The exponent 
— 1 of q is minimal; if it is replaced by any number x < —1, then there exist a £ IR\Q 
such that only finitely many integers q satisfy ||ga|| < q^ . 

Hurwitz 6] proved that the continued fraction algorithm finds, for every a € M\Q, 
an infinite sequence of increasing integers g„ with 

II II 1 -1 
\\Qna\\ <-j^q^ . 

This result is sharp: if the constant is replaced by any smaller one, then the 
statement becomes incorrect. Legendre |12j showed that the continued fraction 
algorithm finds all good approximations, in the sense that if 

II II 1 -1 
mnK-q , 

then q is one of the found by the algorithm. 

As to the generalization of approximations in higher dimensions Dirichlet proved 
the following theorem; See Chapter II of [T6] . 

Theorem 2. Let an n x m matrix A with entries aij € M \ Q be given and suppose 

that 1, Oil, . . . , are linearly independent over Q for some i with 1 < i < n. 

There exist infinitely many coprime m-tuples of integers gi, . . . , such that with 

q = max jgJ > 1, we have 
j 

(3) max \\qi aa H h a„„|| < q^ . 

i 

The exponent of q is minimal. 

Definition 4. Let an n x m matrix A with entries e M \ Q be given. The 
Dirichlet coefficient of an m-tuple qi, . . . ,qm is defined as 

g^maxllgi H ^qmOimW ■ 
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The proof of the theorem does not give an efficient way of finding a series of ap- 
proximations with a Dirichlet coefficient of f . For the case m = 1 the first multi- 
dimensional continued fraction algorithm was given by Jacobi [7]- Many more fol- 
lowed, see for instance Perron Brun |3j, Lagarias pjj and Just [H]. Brentjes [5] 
gives a detailed history and description of such algorithms. Schweiger's book [T7] 
gives a broad overview. For n = 1 there is amongst others the algorithm by Fer- 
guson and Forcade [5]. However, there is no efficient algorithm that guarantees to 
find a series of approximations with Dirichlet coefficient smaller than 1. In 1982 the 
LLL-algorithm for lattice basis reduction was published in 13] . The authors noted 
that their algorithm could be used for finding Diophantine approximations of given 



rationals with Dirichlet coefficient only depending on the dimension; see ( 14 ) 



Just [5] developed an algorithm based on lattice reduction that detects Z-linear 
dependence in the at, in this case m = 1. If no such dependence is found her 
algorithm returns integers q with 

/ n \ 1/2 
max||ga,|| < c ( E q-^^^^"^''-^'>'> , 

where c is a constant depending on n. The exponent — l/(2?i(n— 1)) is larger than 
the Dirichlet exponent — 1/n. 

Lagarias [TU] used the LLL-algorithm in a series of lattices to find good approx- 
imations for the case m = 1. Let ai,...,a„ e Q and let there be a QeN 

with 1 < Q < iV such that maxjlQajH < e. Then Lagarias' algorithm on input 

j 

ai,...,a„ and N finds in polynomial time a q with 1 < q < 2^ N such that 
maxllgajll < V5n2^~e. The main difference with our work is that Lagarias fo- 

cuses on the quality | 1 1 , while we focus on Dirichlet coefficient g " | 1 1 . Besides 
that we also consider the case m > 1. 

The main result of the present paper is an algorithm that by iterating the LLL- 
algorithm gives a series of approximations of given rationals with optimal Dirichlet 
exponent. Where the LLL-algorithm gives one approximation our dynamic al- 
gorithm gives a series of successive approximations. To be more precise: For a 
given n x m-matrix A with entries a^- G Q and a given upper bound gmax the 
algorithm returns a sequence of m-tuples gi , . . . , such that for every Q with 
2 ~ < Q < Qniax One of these m-tuples satisfies 

max \ qj\ < Q and 

j 

(,n + n + 3)(m + n) -,„ 

max||(jiaiiH hgr„a„„||<2 4„ Q" . 

i 

The exponent —m/n of Q can not be improved and therefore we say that these 
approximations have optimal Dirichlet exponent. 

Our algorithm is a multi-dimensional continued fraction algorithm in the sense 
that we work in a lattice basis and that we only interchange basis vectors and add 
integer multiples of basis vectors to another. Our algorithm differs from other multi- 
dimensional continued fraction algorithms in that the lattice is not fixed across the 
iterations. 

In Lemma [25] we show that if there exists an extremely good approximation, our 



algorithm finds a very good one. We derive in Theorem 34 how the output of 



our algorithm gives a lower bound on the quality of possible approximations with 
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coefficients up to a certain Hmit. If the lower bound is positive this proves that 
there do not exist Unear dependencies with all coefficients qi below the limit. 

In Section |4] we show that a slightly modified version of our algorithm runs in 
polynomial time. In Section [5] we present some numerical data. 

2. Lattice reduction and the LLL-algorithm 

In this section we give the definitions and results that we need for our algorithm. 

Let r be a positive integer. A subset L of the r-dimensional real vector space W is 
called a lattice if there exists a basis bi, . . . ,br of E'' such that 

r f r \ 

L = ^ Zfe, = i ^ zA; z^eZ{l<i<r)^ . 
1=1 1 1=1 J 

We say that 61, . . . , 6^ is a basis for L. The determinant of the lattice L is defined 
by I det(&i, . . . , br)\ and we denote it as det(L). 

For any linearly independent bi, . . . ,br € M'' the Gram-Schmidt process yields an 
orthogonal basis 6^, . . . , 5* for W, by defining inductively 

i-l 

(5) b* =bi~ ^ ^Ji,ijb* for 1 < i < r and 

where ( , ) denotes the ordinary inner product on W . 
We call a basis bi, . . . ,br for a lattice L reduced if 

l^iil < ^ for 1 < j < i < r 

and 

\b*+^,u-lbUf <l\bUf forl<z<r, 
where denotes the Euclidean length of x. 
The following two propositions were proven in |14j . 

Proposition 6. Let bi, . . . ,br be a reduced basis for a lattice L inMJ . Then we 
have 

{^) |6l|<2(-l)/4(det(L))^/^ 

(a) |6i P < 2'-^ for every x eL,x^O, 

r 

{lii) < 2''('--i)/4det(L). 

1=1 

Proposition 7. Let L (Z XT be a lattice with a basis 61, b2, ■ ■ ■ , br, and let F d M., 

F > 2, be such that \bi\'^ < F for 1 < i < r. Then the number of arithmetic 
operations needed by the LLL-algorithm is 0{r'^ log F) and the integers on which 
these operations are performed each have binary length 0{rlogF). 



In the following Lemma the approach suggested in the original LLL-papcr for find- 
ing (simultaneous) Diophantine approximations is generalized to the case m > 1. 
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Lemma 8. Let an nx m-matrix A with entries aij in M and an e G (0, 1) be given. 
Applying the LLL-algorithm to the basis formed by the columns of the (m + n) x 
(m + n) -matrix 



(9) 



B = 



1 

1 

... 

... 



an 
a2i 



a2r, 



1 




c 







with c = ^2 ^4 yields an m-tuple 51, . . . , (7„i G Q with 

(10) max|gj| < 2~ ^~ £~ and 

i 



(11) 



maxllgiOiiH \- qmaim\\ < £■ 



Proof. The LLL-algorithm finds a reduced basis 61, ... , bm+n for this lattice. Each 
vector in this basis can be written as 

giaii H h qmaim - Pi 



QlClnl + • ■ ■ + qmO-nm — Pn 

cqi 
cqni 

with Pi G Z for 1 < i < n and qj G Z for I < j < m. 

Proposition [gJ^I) gives an upper bound for the length of the first basis vector, 

|5i| < 2'''t^c^. 
From this vector bi we find integers qi , . . . , , such that 
(12) 
(13) 



max llgiOii 



max|(7j| < 2 ^4 0™+" and 
i 

< 2 — 3 — C^. 



qma 



Substituting c = ^2 "^4' ^ " gives the results. 



□ 



From equations ( 12 ) and ( 13 ) it easily follows that the m-tuple qi, . . . ,q„i satisfies 
(14) 



max||giaii-| hgr„a^„||<2 q" 



where g = max|gj|, so the approximation has a Dirichlet coefficient of at most 



, (m + ,^-l)(r7l + ,^) 
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3. The Iterated LLL-algorithm 

We iterate the LLL-algorithm over a series of lattices to find a sequence of approx- 
imations. We start with a lattice determined by a basis of the form After the 
LLL-algorithm finds a reduced basis for this lattice, we decrease the constant c by 
dividing the last m rows of the matrix by a constant greater than 1. By doing so, 
e is divided by this constant to the power We repeat this process until the 



upper bound (10 1 for q guaranteed by the LLL-algorithm exceeds a given upper 



bound Qmax- Motivated by the independence on e of ( 14) we ease notation by fixing 



Define 

(15) k' := 



(m -|- n — l)(r7i -|- n) TOlog2(?n 



4n 



Iterated LLL-algorithm (ILLL) 

Input 

An n X 

An upper bound <7max > 1- 



An n X m-matrix A with entries a^- in . 



Output 

For each integer k with 1 < fc < fc', see (15), we obtain a vector q{k) E Z™ with 



(16) max|g,(fc)| < 2 ""^"'4^'""^'" 2^, 

j 

(17) max||9i(fc)a.a H h g,„(/s) a„„|| < ^. 

Description of the algorithm 

(1) Construct the basis matrix B as given in ^ from A. 

(2) Apply the LLL-algorithm to B. 

(3) Deduce gi, . . . from the first vector in the reduced basis returned by 
the LLL-algorithm. 

(4) Divide the last m rows of B by 2^s^ 



(5) Stop if the upper bound for q guaranteed by the algorithm ( 16 ) exceeds 
gmax; else go to step 2. 



Remark 18. The number 2 ™ in step 4 may be replaced by d m for any real 
number d > 1. When we additionally set e = ^ this yields that 

(19) max|qj-(fc)| < 2 '"'+"-4'2'"'^"' rf^ and 

j 

(20) max||gi(/c)a,a H \- qm.{k)aim\\ < 



-k 



In the theoretical part of this paper we always take d ~ 2 corresponding to our 
choice e = \- 

Lemma 21. Let an n x m-matrix A with entries aij in M and an upper hound 
(Zmax > 1 &e given. The number of times the ILLL-algorithm applies the LLL- 
algorithm on this input equals k' from (15). 
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Proof. One easily derives the number of times we iterate by solving k from the 



stopping criterion ( 16 1 



(rn + n^-l)C21±Ili„fc,J 

< 2 2—, 



□ 



We define 

c{k) = c{k~l)/2— for k>l, where c(l) = c as given in Lemma [8j 

In iteration k we are working in the lattice defined by the basis in ([9| with c replaced 
by c(fc). 



Lemma 22. The k-th output, q{k), of the ILLL-algorithm satisfies (16) and (17) 



Proof. In step k we use c{k) — {2 fe+i^ " _ Substituting c{k) for c in 

equations (12 1 and (13) yields (Il6|) and (17), respectively. □ 



The following theorem gives the main result mentioned in the introduction. The 
algorithm returns a sequence of approximations with all coefficients smaller than 
Q, optimal Dirichlet exponent and Dirichlct coefficient only depending on the di- 
mensions m and n . 

Theorem 23. Let an nx m-matrix A with entries aij in M, and Qmax > I be given. 
The ILLL-algorithm finds a sequence of m-tuples gi, . . . , q„i such that for every Q 
with 2 rsi ^ < Q < (Zmax one of these m-tuples satisfies 

max \ qj\ < Q and 
j 

maxllqifliiH ^q,nairn\\<2 Q - . 



Proof Take fc € N such that 

(24) 2 •2^^<Q<2 ^ -2- 



From Lemma 22 we know that q(k) satisfies the inequality 

(m + ,^ + 3)(,Tl + ,^) . (fc-1) 



max|qj(fc)| < 2 ^ 2^^ < Q. 

j 



From the right hand side of inequality (24) if follows that i < 2* ^ ^4,^* ^ ^ 



From Lemma 22 and this inequality we 



Herive that 



II /1\ /I \ II 1 (Tn + n+3)(m + n) 

max||(7i(fcjaji H \- qm.(k) OimW < ^ < 2 -i- Q " . 

□ 

Proposition |6]^ii) guarantees that if there exists an extremely short vector in the 
lattice, then the LLL-algorithm finds a rather short lattice vector. We extend this 
result to the realm of successive approximations. In the next lemma we show that 
for every very good approximation, the ILLL-algorithm finds a rather good one not 
too far away from it. 
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Lemma 25. Let an nx m-matrix A with entries aij in M, a real number < 6 < 1 
and an integer s > 1 be given. If there exists an m-tuple Si, . . . , s„i with 



(26) s = max Is,- 1 > 2 

] \ m 

and 

(27) max||siaa H hs,„ai„i|| < Ss^ , 

i 

then applying the ILLL- algorithm with 

(28) Qmax > 2 (^ — j S 

yields an m-tuple qi, . . . ,qm with 

(29) max|g,| < 2 ( — ) 

and 

(30) max llqiaji H \- qmai,n\\ < 2^~^/nSs~. 

i 

Proof. Let 1 < fc < fc' be an integer. Proposition |6jii) gives that for each q{k) 
found by the algorithm 

n m 

\\qi{k)aa H hq„i (fc)a^mf +c(fc)2^g,(fc)2 

1=1 3=1 



S 



1 j=i 



From this and ^261 and (271 it follows that 



(31) max||gi(fc)a,a + • • • + 9m(fc)a„„|p < 2"+""^ (^nc^^s^ + c(fc)2ms2) 

Take the smallest positive integer K such that 

(32) c{K) < .f^Ss-'^. 

V m 



We find for step K from (31 1 and (32) 



ma.x\\qi{K)a.a ^ \- qrn{K)aim\\ < 2 2 y^Ss 



which gives ( 30 ) . 



We show that under assumption (28) the ILLL-algorithm makes at least K steps. 
We may assume K > 1, since the ILLL-algorithm always makes at least 1 step. 
From Lemma [2l] we find that if <7max satisfies 

A'n (m + n-l)(m + r.) 

(/max ^ ^ ^ ^"^ , 

then the ILLL-algorithm makes at least K steps. Our choice of K implies 

, s _ c(l) _ 2 4,„ m+2L 

cl-"- ^) - ^ (,„+„)(fc^2) - ^ (,„+,>)(g-2) \j rn " ' 



and we obtain 



2— < 2' 



}2t ( rn \ 
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From this we find that 

9 , , . . . . y , H 



m^+m(n-l) + 4n / TO \ 



S 



is a satisfying condition to guarantee that the algorithm makes at least K steps. 

Furthermore, either 2 ^ ' ^ < or X = 1. In the former case we 

find from ^ that 

ma.x\q,{K)\<2'^c{K)^^ < 2^ ( ^)^^^ s . 

j \no^ J 

In the latter case we obtain from (TTll 



max|9j(l)| < 2™^" = 2 ' 2 '""'^'J,!^'" 



and, by (26), 



2^^2 — ^ — =2^i^2"2 — — <2^i^2™ ^ 
We conclude that for all if > 1 



max|g,(if)|<2 - ( — j 



□ 



Note that from (29) and (30) it follows that 



(33) q"^ ia&yi\\qiaii + ■ ■ ■ + q^a^mW < 2™ +"<"",^'+""+"" ^2(r+„) (^52^ 2(™+„) ^ 

where again g = max 
j 

Theorem 34. Let an n x m-matrix A with entries aij in R and gmax > 1 be 
given. Assume that 7 is such that for every m-tuple qi,...,qm returned by the 
ILLL-algorithm 

(35) q~ max Hgiflii + . . . qm0.im\\ > 7, where q = max \qj\. 

i j 

Then every m-tuple Si, . . . , with s = maxj \ sj\ and 



(„ + „-!). 



2 4?;^ — <s < 2 — q 

m V TO / 



max 



maxllsiOii H HSmflim II > 

i 

with 

^ -(m + ,i)(„i^+r,i(3r.-l) + 4,» + 2,i^) _„ _1 rr^±n 

(36) (5 = 2 4„2 TO^n~7~ 



Proof. Assume that every vector returned by our algorithm satisfies ( 35 ) and that 
there exists an TO-tuple si, . . . , Sm with s ~ max^ |sj| such that 

/n(5^\ ^<^^ ,„^+^(„-i)+4„ 

2 4„ <S<2 4™ Qmax 

\ m J \ m J 

and s max ||sia,;i + • • • + SmaimW < <5. 
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From the upper bound on s it follows that Qmax satisfies (28). We apply Lemma 25 
and find that the algorithm finds an m-tuple qi, . . . ,qm that satisfies ( 33 ) . Substi- 



tuting S as given in ( 36 ) gives 



g" max||(jiaii 



< 



7, 



which is a contradiction with our assumption. 



□ 



4. A POLYNOMIAL TIME VERSION OF THE ILLL- ALGORITHM 

We have used real numbers in our theoretical results, but in a practical implemen- 
tation of the algorithm we only use rational numbers. Without loss of generality we 
may assume that these numbers are in the interval [0, 1]. In this section we describe 
the necessary changes to the algorithm and we show that this modified version of 
the algorithm runs in polynomial time. 

As input for the rational algorithm we take 



• the dimensions m and n, 

• a rational number e € (0, 1), 



an integer M that is large compared to ^" 



m+n 



loge, 



• an n X m-matrix A with entries < a^^ < 1, where each Uij — ^ for some 



integer pij, 
• an integer q„ 



< 2 



M 



When we construct the matrix B in step 1 of the ILLL-algorithm we approximate 
c as given in (|9]) by a rational 



■)M 



2M 



1 

2M ■ 



(37) 

Hence c < c < c 
In iteration k we use a rational c{k) that for fc > 2 is given by 

2*-^c(fc- 1)2-"^ 



c(fc) 



and c(l) = c as in (37 1 



and we change step 4 of the ILLL-algorithm to 'multiply the last m rows of B by 
c{k — l)/c(fc)'. The other steps of the rational iterated algorithm are as described 
in Section ini 



4.1. The running time of the rational algorithm. 

Theorem 38. Let the input he given as described above. Then the number of 
arithmetic operations needed by the ILLL-algorithm and the binary length of the 
integers on which these operations are performed are both bounded by a polynomial 
in m, n and M . 



Proof. The number of times we apply the LLL-algorithm is not changed by ratio- 

)er of steps fc' from Lemma 21 

n — 1) (to + n) TO log2 (/max 1 \ TuM 



nalizing c, so we find the number of steps fc' from Lemma 21 
(to 



fc' 



An 



< 
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It is obvious that steps 1, 3, 4 and 5 of the algorithm are polynomial in the size 
of the input and we focus on the LLL-step. We determine an upper bound for the 
length of a basis vector used at the beginning of an iteration in the ILLL-algorithm. 

In the first application of the LLL-algorithm the length of the initial basis vectors 
as given in (|9| is bounded by 

l&ip < max {l, afj + • • • + af^j + mc^} < m + n, for 1 < i < m + n . 
where we use that < a^j < 1 and c < 1. 

The input of each following application of the LLL-algorithm is derived from the 
reduced basis found in the previous iteration by making some of the entries strictly 
smaller. Part (ii) of Proposition |6] yields that for every vector bi in a reduced basis 
it holds that 

|b.P<2 """""""""^' (det(£))^ n l^r'- 

j = l J't^* 

The determinant of our starting lattice is given by c™ and the determinants of all 
subsequent lattices are strictly smaller. Every vector bi in the lattice is at least 
as long as the shortest non-zero vector in the lattice. Thus for each i we have 
l^iP > Combining this yields 

,9 (m + ,i + 2AJ)(m + ,.-l) , 

l^il < 2 2 c^"' < 2" 

for every vector used as input for the LLL-stcp after the first iteration. 
So we have 



(39) < max {m n , 2 """""^"""-""'^' } 



2 2 



for any basis vector that is used as input for an LLL-step in the ILLL-algorithm. 

Proposition [t] shows that for a given basis bi, . . . ,bm+n for Z'"+" with € M, 
F > 2 such that < F for 1 < i < m + n the number of arithmetic operations 
needed to find a reduced basis from this input is 0{{m + n)^\ogF). For matrices 
with entries in Q we need to clear denominators before applying this proposition. 
Thus for a basis with basis vectors < F and rational entries that can all be 
written as fractions with denominator 2^^ the number of arithmetic operations is 
0((m-hn)4log(22*^F)). 

Combining this with (|39|) and the number of steps yields the proposition. □ 



4.2. Approximation results from the rational algorithm. Assume that the 
input matrix A (with entries a.y — |^ € Q) is an approximation of an n x m-matrix 

A (with entries atj S K), found by putting aij — - — p?^- In this subsection we 
derive the approximation results guaranteed by the rational iterated algorithm for 
the e M. 
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According to ( 12 ) and ( 13 ) the LLL-algorithm applied with c instead of c guarantees 



to find an m-tuplc qi, . . . , q„i such that 



max Wqian 



q — max I g, | 
i 

and 



„(i:i±izil)(zzi±ii) - 

< 2 £ 



II < 2-^(^(2--^e) 

(m + n-l)(m + n)-4A/ni 
< e + 2 4(™ + „) 



1 ^ '"+" 
2^ 



the last inequality follows from the fact that {x + y)" < a;" + if a < 1 and 
x,y >0. 



For the aij we find that 

max llgia^i H h qmaimW < max ||giaa H h gma™|| + 77192 



m + n-l Mm -n {m + n - 1) jm + r, 

<e + 2 '"+" + 777£ " 2 'I™ 



On page [9] we have chosen M large enough to guarantee that the error introduced 
by rationalizing the entries is negligible. 

We show that in every step the difference between c(/c) and c(fc) is bounded by w • 



Lemma 40. For each integer k > 0, 



m < m < c{k) + ^ ^ 2-^ < c{k) 



i=0 



2 

2^' 



Proof. We use induction. For k — we have c(0) = - — — and trivially 



c(0) < c(0) < c(0) 



1 

w 

fe-1 



Assume that c(fc — 1) < c(fc — 1) < c(fc — 1) + — ^ ^ 2 ' ™ ' and consider c{k) 



i=0 



From the definition of c(fc) and the induction assumption it follows that 



c(fc- 1) 2-^=^2*^ 



and 



c{k - 1) 2- 



< 



< 



c(fc- 1) 



2Af 

1 v^fc-1 



:=(fc-l)+2^Eto2 



2^ 
fe 



+ 



)M 



i=0 



Finally note that 2 ' ™ ' < 2 for all fc. 



i=0 



□ 
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One can derive analogues of Theorem[23] Lemma[25]and Theorem 34 for the polyno- 
mial version of the ILLL-algorithm by carefully adjusting for the introduced error. 
We do not give the details, since in practice this error is negligible. 



5. Experimental data 

In this section we present some experimental data from the rational ILLL-algorithm. 
In our experiments we choose the dimensions m and n and iteration speed d. We fill 
the mxn matrix A with random numbers in the interval [0, 1] and repeat the entire 
ILLL-algorithm for a large number of these random matrices to find our results. 
First we look at the distribution of the approximation quality. Then we look at the 
growth of the denominators q found by the algorithm. 



5.1. The distribution of the approximation qualities. For one-dimensional 
continued fractions the approximation coefficients Ok are defined as 

r, 2 Pk 

Qk 

where Pk/qk is the nth convergent of a. 

For the multi-dimensional case we define Qk in a similar way 
(41) Gfc = g(fc)^ max||(7i(A;)a.a H hg,„(fc)ai„ 



The one- dimensional case m = n = 1. In [T| it was shown that for optimal contin- 
ued fractions for almost all a one has that 



lim 



N- 



i# {1 < n < TV : e„(x) <z}^ F{z), where 



logG' 



F{z) = I VT^ + log{G'-^^q^) 
logG 



75 



< z < 



1, 



< z < 1, 



where G= 



As the name suggests, the optimal continued fraction algorithm gives the optimal 
approximation results. The denominators it finds, grow with maximal rate and all 
approximations with 6 < ^ are found. 

We plot the distribution of the 0's found by the ILLL-algorithm for m = n = 1 
and d = 2 in Figure [T] The ILLL-algorithm might find the same approximation 
more than once. We see in Figure [l] that for d = 2 the distribution function 
differs depending on whether we leave in the duplicates or sort them out. With the 
duplicate approximations removed the distribution of Q strongly resembles F{z) of 
the optimal continued fraction. The duplicates that the ILLL-algorithm finds are 
usually good approximations: if they are much better than necessary they will also 
be an admissible solution in the next few iterations. 



For larger d we do not find so many duplicates, because the quality has to improve 
much more in every step; also see Figure [2] for an example with d = 64. 



DIOPHANTINE APPROXIMATIONS WITH PRESCRIBED QUALITY 



// 



Optimal CF 



d=2 with duplicates 



Figure 1. The distribution function for from ILLL with 
m = n = 1 and d = 2, with and without the duphcate approxi- 
mations, compared to the distribution function of 6 for optimal 
continued fractions. 




OptiinalCF 



d=64 duplicates 



Figure 2. The distribution function for @ from ILLL with 
m = n = 1 and d = 64, with and without the duphcate approx- 
imations, compared to the distribution function of O for optimal 
continued fractions. 



now on we remove duplicates from our results. 




Figure 3. The distribution function for Q from ILLL (with du- 
plicates removed) with m = n = 1 and various values of d. 
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5.2. The multi-dimensional case. In this section we show some results for the 
distribution of the 0's found by the ILLL-algorithm. For fixed m and n there also 
appears to be a limit distribution for as d grows. See Figure |4] for an example 
with TO = 3 and n = 2, and compare this with Figure [3j In this section we fix 
d = 512. 



m = 3 n = 2 




0.0 0.2 04 06 OB i.O 



Figure 4. The distribution function for 9 from ILLL with to = 3 
and n = 2 for d = 2, 8, 128 and 512. 

In Figure [5] we show some distributions for cases where either to or n is 1. 




Figure 5. The distribution for Q from ILLL when either to = 1 
or n = 1. 

In Figure |6] we show some distributions for cases where m = n. 

Remark 42. Very rarely the ILLL-algorithm returns an approximation with > 1, 
but this is not visible in the images in this section. 

5.3. The denominators q. For regular continued fractions, the denominators 
grow exponentially fast, to be more precise, for almost all x we have that 

lim ql^'' = ei2Tog2 ^ 

see Section 3.5 of [1]. 

2 , 2 

For nearest integer continued fractions the constant y^i^^2 ^® replaced by ]^2k)gG 
with G = ^.^^ ■ For multi-dimensional continued fraction algorithms little is known 




Figure 6. The distribution of 6 from ILLL when m — n. 

about the distribution of the denominators qj . Lagarias defined in ^ the notion of 
a best simuhaneous Diophantine approximation and showed that for the ordered 
denominators 1 = gi < (72 < • • • of best approximations for ai, . . . , a„ it holds that 

hm inf q\!^ > 1 + — ^--r- 



m=l,n=l,d = 2 
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A. 



Wlrrr>rh 



12 14 16 18 20 



12 14 16 18 20 



Figure 7. Histograms of e for various values of m,n and 

d. In these experiments we used qmax = 10'*" and repeated the 



ILLL-algorithm [^^t^J times, with k' from Lemma 



21 



We look at the growth of the denominators q — maxj \qj\ that arc found by the 
ILLL-algorithm. Dirichlct's Theorem [2] suggests that if q grows exponentially with 
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a rate of m/n, then infinitely many approximation with Dirichlet coefficient smaller 
than f can be found. In the iterated LLL-algorithm it is guaranteed by ( 16 ) that 
q{k) is smaller than a constant times . Our experiments indicate that q{k) 

kn . ^ log <?fc . . I I 

is about , or equivalently that e is about d; see Figure M which gives a 

m log 

histogram of solutions that satisfy e = x. 
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